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We study the tunneling through an arbitrary number of finite rectangular opaque barriers and 
generalize earlier results by showing that the total tunneling phase time depends neither on the 
barrier thickness nor on the inter-barrier separation. We also predict two novel peculiar features of 
the system considered, namely the independence of the transit time (for non resonant tunneling) and 
the resonant frequency on the number of barriers crossed, which can be directly tested in photonic 
experiments. A thorough analysis of the role played by inter-barrier multiple reflections and a 
physical interpretation of the results obtained is reported, showing that multibarrier tunneling is a 
highly non-local phenomenon. 
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I. INTRODUCTION 

A renewed interest in a typical quantum phenomenon 
such as the tunnel effect has been recently achieved due to 
a long series of experiments aimed to measure the tunnel- 
ing transit time (for reviews see, for istance, [1]). While 
such experiments involving electrons are usually difficult 
to realize (mainly due to the smallness of the electron 
de Broglie wavelength at usual temperatures) and even 
of uncertain interpretation, the observations on photonic 
tunneling [2]- [6] has by now provided clear data on this 
subject. Despite the different phenomena studied in sev- 
eral experiments (undersized waveguides, photonic band 
gap, total internal reflection) and the different frequency 
ranges for the ligth used (from the optical to the mi- 
crowave region), all such experiments have shown that, 
in the limit of opaque barriers, the transit time to travel 
across a barrier of width a is usually shorter than the cor- 
responding one required for real (not evanescent) propa- 
gation through the same region of width a. This result 
can be interpreted [1] in terms of a superluminal group 
velocity i> gr > c which, however, does not violate Einstein 
causality, since the signal velocity relevant for that [7] is 
never measured. Nevertheless we prefer to look at the 
experimental result as an observation of the simple Hart- 
man effect [8]: for opaque barriers the tunneling phase 
time is independent of the barrier width. Although sev- 
eral definitions of the tunneling time (also related to the 
different experimental setups used) exist [1] and a gen- 
eral consensus on this is still lacking, it seems that all 
the experimental results can be successfully interpreted 
in terms of phase time [9]. 

Further light has been put on the problem by recent ex- 
periments involving double barrier penetration [10]. In 
fact, while the above effect has been confirmed in such 



a system too (far from the resonances of the structure), 
observations show that the transit time is also indepen- 
dent of the separation distance between the barriers (sup- 
posed to be thick). This peculair phenomenon has been 
studied theoretically in Ref. [11], where the authors have 
provided a straigthforward generalization of the Hartman 
effect for double barrier tunneling. 

Convinving qualitative explanations of these two findings 
(namely that the tunneling phase time is independent of 
the barrier thickness as well as of the inter-barrier sep- 
aration for opaque barriers) have been reported. When 
considering a given wavepacket entering into a potential 
barrier region, a reshaping phenomenon occurs in which 
the travelling edge of the pulse is preferentially attenu- 
ated that the leading one, thus simulating a group veloc- 
ity greater than c [1] , [3] . In practice the Hartman effect 
in the tunneling through a thick barrier is explained from 
the fact that under the barrier no phase accumulates, 
and the entire phase shift comes only from the bound- 
aries, thus being substantially independent of the thick- 
ness [12]. Furthermore, when two barriers are present, 
the transit time independence on the barrier separation 
can, instead, be understood in terms of an effective accel- 
eration of the forward travelling waves in the inter-barrier 
region, which arises from a destructive interference be- 
tween the two barriers [11]. 

Further noticeable results have been recently achieved 
in Ref. [12], where it has been shown that a wavepacket 
travels in zero time a region with N arbitrary 5— function 
barriers. 

In this paper we extend all these findings by consider- 
ing the case of N successive opaque barriers with finite 
widths and heigths. While we confirm all previous re- 
sults, we generalize them by showing that some pecu- 
liar tunneling properties are independent of the number 
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FIG. 1. Potential barrier V(x) with TV equally spaced rect- 
angular bumps of given heigth and width. 



following form *: 



ip(x) = < 



f ip2i(x) (i-l)L<x<(i-l)L + a 
(< = 1,2 JV) 



i>2i+l(x) 



otherwise 
(i = 0,l,2,...,TV) 



(2) 



with: 



^(x) = e lkx + Re 



-ikx 



(3) 



of the barriers crossed (Sect. II). Furthermore, in or- 
der to establish a quantitative interpretation of the in- 
volved phenomena, in Sect. Ill we study the role of mul- 
tiple reflections in double barrier tunneling and show how 
strongly the total tunneling phase time depends on non- 
local effects. Finally in Sect. IV we discuss the results 
obtained and give our conclusions. 

In view of the formal analogu [13] between the 
Schrodinger equation and the electromagnetic Hclmholtz 
equation, our study applies to matter particle tunnel- 
ing as well as to evanescent propagation of photonic 
wavepackets. This is a straigthforward cosequence of the 
fact that in both cases the starting point is basically the 
same (in our case it is Eq. (2)) [9], on interchanging 
the roles of angular frequency lo and wavevector k into 
the corresponding ones of energy E and momentum p 
through the Planck - de Broglie relations. Thus, through- 
out this paper, we use indifferently particle or wave termi- 
nology unless the meaning of what we are doing becomes 
unclear. 



II. TUNNELING THROUGH TV SUCCESSIVE 
BARRIERS 

Let us consider a wavepacket moving along the x— axis 
and entering at x — into a region with a potential 
barrier V(x) as depicted in Fig. 1: 



V(x) 



V (i-l)L<x<(i-l)L + a 







otherwise 



(1) 



( i=l,2,. . . ,N). For the sake of simplicity we choose the 
heigth Vo of the potential barriers, as well as the width a 
of each barrier, to be the same for all TV rectangular bar- 
riers. We further assume equally spaced barriers, L — a 
being the inter-barrier distance. 

The propagation of the wavepacket through the bar- 
riers is described by a scalar field ip representing the 
Schrodinger wavefunction in the particle case or some 
scalar component of the electric or magnetic field in the 
photonic case. This is the solution of the Schrodinger 
equation or the Hclmholtz equation with potential or re- 
fractive index in Eq.(l) and, in both cases, it takes the 



(i = 1,2 JV) (4) 

(i = l,2 JV-1) (5) 



i; 2N+1 (x)=Te lk ^ N - 1 ^ 



(6) 



As noted in Ref. [9], the explicit dependence on the fre- 
quency of the (real) wavevector k in the barrier-free re- 
gions and imaginary wavevector i\ in the barrier ones 
enters only in the final expression for the phase time. As 
long as possible we do not use a particular dispersion 
relation in order to draw general features which are com- 
mon to the particle and to the wave case. 
The 4TV unknown coefficients R,T, Ai, Bi are obtained 
from the AN matching conditions for the function ip and 
its derivative ip' at the discontinuity points x + (i — 1)L, 
x = (i — 1)L + a of the potential: 



(7) 



(i = 1,2,..., N). Note that the quantities R and T have 
the meaning of (total) reflection and transmission coeffi- 
cient from the TV— barrier system, respectively, and sat- 
isfy the unitarity condition: 



^2i-l {X = 


(i-l)L) 


= i>2i (x = 


(i 


-1)L) 


^2i-\ ( x = 


(i-l)L) 




(i 


-l)L) 


1p2i (x = (i - 


l)L + a) 


= i>2i+i (x 




(*-l)L + a) 


i>2t (x = {i- 


l)L + a) 


= i>2i+l(x 




(*-l)L + o) 



\R\ 2 + \T\ 



(8) 



* Obviously, the physical field is represented by a wavepacket 
with a given spectrum in uj: 



dio rj(ui) ip(x) e 



where 77(01) is the envelope function. Keeping this in mind, for 
the sake of simplicity we deal with only stationary solutions 
as in Eq. (2). 
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We have produced a Mathematica symbolic code in order 
to obtain explicit analytic expressions for all the coeffi- 
cients appearing in Eqs. (3)-(6). However, we here report 
only the interesting result obtained for the transmission 
coefficient T(N) for an N— barrier system in the opaque 
barrier approximation %o >> 1. In this limit the quan- 
tity T(N) can be factorized in the following way: 



T(N) 



ika 



Co • £{N) ■ F{N) 



(9) 



Co 
£(N) 



4%xk 
(k + ix) 2 

r -^„iN 



2 x k 



2xfccos k(L — a) — (k 2 — x 2 ) sinfc(L — a) 



N-l 



Note that only the real terms £ and T depend on a, L, N, 
while the complex factor Co does not. As a consequence, 
since the tunneling phase time t is defined as: 



du 



(10) 



and the quantity 



arg{T(A0e 4fca } - arg 



Ai X k 
(k + i X ) 2 



arctan 



2 X k 



(11) 



is independent of a, L, N, we arrive at the general con- 
clusion that the tunneling phase time for a system of N 
opaque barriers depends neither on the barrier width and 
inter-barrier distance nor on the number of the barriers. 
Let us now discuss the effects of the real terms in Eq.(9) 
on the tunneling probability P T (N) = \T(N)\ 2 : 



Pt(N) 



4%fc 
k 2 + x 2 



2 x k 



i 2N 



2xk cos k(L — a) — (k 2 — X 2 ) sin k(L — a) 



2(N-1) 



(12) 



We easily recognize that the last factor in Eq.(12), com- 
ing from the term T{N) is responsible of the resonance 
structure of the transmission probability. The factor 
F(N) is, of course, absent in the case of only one barrier, 
i.e. N = 1 or N ^ 1 but L = a. However, no reso- 
nance can occur even in the particular case in which the 
inter-barrier distance is tuned in a way that: 



L — a = 



~~k 



(i/ = 0,1,2,...) 



(13) 



In this case, waves moving forward and backward in the 
inter-barrier regions interfere between them such that no 
resonance takes place. 

The resonance condition for the tunneling probability is, 
from Eq. (12), the following: 



tanfc(i — a) 



%xk 

k 2 — x 2 



(14) 



It is wortwhile to observe that Eq.(14) does not depend 
on N, so that the resonant frequency is the same irre- 
spective of the number of barriers to be crossed. 
Finally we point out an intriguing consequence of the 
resonance condition. Let us write Eq.(14) as follows: 



tan</> • t&nk(L — a) 



(15) 



where <f> is given in Eq.(ll), and take the derivative of 
Eq. (15) with respect to the angular frequency u>. By 
using Eq. (10) we easily arrive at the following relation: 



T + T 







(16) 



where r is the (phase) time for travelling the inter- 
barrier distance L — a in vacuum. Keeping in mind that 
the total tunneling time has the same value of the tun- 
neling time for crossing only one barrier (see above), we 
see that, when resonant tunneling takes place, the total 
time required to cover the distance L (one barrier length 
a plus one inter-barrier distance L — a) is zero. 



III. MULTIPLE REFLECTIONS AND 
NON-LOCALITY 

In order to have a physical interpretation of the results 
obtained previously, we now consider the effect of single 
barriers on the propagation of the wavepacket through 
the entire N— barrier system, by invoking the superposi- 
tion principle. For the sake of simplicity, we will study 
the case of a system of two opaque barriers. 



A. Partial coefficients 

For N = 2, in the barrier- free regions, Eqs. (3)- (6) reduce 
to the following: 



-ikx 



fa(x) = e l * x + Re 
fo(x) = A 3 e ikx +B 3 e 
iP b {x)=Te lk{x - L) 



-ikx 



(17) 



where the explicit expressions for the coefficients are re- 
ported in Appendix A. Let us now denote with R\,T\ 
and i?2,T 2 the (partial) reflection and transmission co- 
efficients of the first and second barrier, respectively. In 
the region with x < the reflected wave is described by 
the term: 



Re- lkx = R ie - lkx + B 3 T ie - ikx , (18) 
while for x > L + a the transmitted one is described by: 
Te ik(x-L) = AzT 2 e ik{x ~ L) . (19) 
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By taking into account multiple reflections from the two 
barriers in the region with a < x < L, we see that the 
forward travelling wave is described by the term: 



A 3 e lkx = T 1 \l + R 1 R 2 + (R 1 R 2 f + . 

while the backward one is described by: 

B 3 e~ ikx = A 3 R 2 e" ik ^ L ^ 
Then, by introducing the quantity: 

{RiR2y = — 1 



'ikx 



s = J2 ( R i R i) 1 



1 — R\R 2 



(20) 



(21) 



(22) 



which accounts for multiple reflections, from Eqs.(18)- 
(21) we obtain: 



R — R\ + B 3 Ti 

T = A 3 T 2 

A 3 = TiS 

B 3 = A 3 R 2 e ikL 



(23) 



By solving these equation with respect to the partial re- 
flection and transmission coefficients, we get: 



Ri 



R 2 = 



T 2 = — 



R- 


A 3 B 3 


1 - 


-Bl 


A 3 - 


- B 3 R 


1 - 




B 3 


—ikL 


T 3 e 




T 




M 





(24) 



In the opaque barrier limit x a >> 1j f° r the second bar- 
rier we obtain: 



R2 

T 2 



Rob e tkL 
Tqb e lkL 



while for the first barrier: 



Ri — Rob + Rq + Rr 
T\ = Tqb +Tq + T r 



(25) 



(26) 



where: 



^OB 



Toe 



k + i\ 

Ai X k 
(k + ix)' 



Aixk 
(k + ix) 2 ' 



-2 X a 



(27) 



-e~ ika e~ xa 



are the reflection and transmission coefficients corre- 
sponding to a one-barrier system (TV = 1) and: 



Rq 



k-ix 
k + ix 



■p2 e 2ik{L-a) e -2 X a 



R R = [ T 2 e lkL e- 2 * a 



k + ix 

k — ix 
k + ix 

k-ix 
k + ix 



(28) 



jT e 2ik(L-a) e -ikL e -xa 



For future reference, we also consider the partial coeffi- 
cients Ri, Tj 3 , i?2i T° in the approximation of no multiple 
reflections, as considered in Ref. [11] t. These are ob- 
tained from Eqs.(23) by setting S = 1. We have 



Ri — Rob + Rq 
Ti = T b + Tq , 

while R 2 ,T^ are the same as in Eqs.(25). 



(29) 



B. Unitarity conditions 

The interpretation of the quantities R\,T\ and R 2 , T 2 as 
reflection and transmission coefficients of the first and 
second barrier is derived from the unitarity conditions 
satisfied by these coefficients. In fact, since |i?| 2 + |T| 2 = 
1 and I^obI 2 + |7ob| 2 = 1, wc find that: 



|i?i| 2 + |T!| 2 = 1 

r 2 | 



Ii? 2 | 2 + IT 2 | 2 = 1 



(30) 



It is easily recognizable as well that, assuming no multiple 
reflection, the total probability for scattering from the 
first barrier is lower than 1 5 : 



R\\ 2 + \T?\ Z = l-f z e 



i0|2 



(31) 



^The authors of Ref. [11] have considered the case of no 
multiple reflections and, moreover, they also neglected the 
second term B3T1 in the first equation in (23) corresponding 
to backward waves in the x < region transmitted from the 
first barrier, reflected from the second one and again transmit- 
ted from the first barrier. While their parametrization of the 
wavefunction is, of course, permitted and leads to correct re- 
sults, nevertheless the partial coefficients they obtained have 
no direct physical meaning, as we will show below. 

*In the approximation considered in Ref. [11] the quantity 
Ri in Eq.(29) should be replaced by the following one: 



Ri = -Rob + 



k — i\ 4:ixk 



k + ix (k + ix) 2 



ik(L-a) g-2x<2 



§ Instead, by using the parametrization of Ref. [11] , we obtain 
an unphysical scattering probability greater than 1, 
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this revealing that something has been forgotten. Obvi- 
ously, multiple reflections are the missing term and it 
is worth to observe that the probability for this phe- 
nomenon to occur, which from Eq.(31) we deduce to be 
.F 2 e -2xa , is given by: 



\Rr\' 



-p2 e -2 X a 



(32) 



Thus the quantities i?R and Tr, that must be added to 
the no multiple reflection coefficients i?° and T® in order 
to obtain the complete ones i?i and T\ respectively, can 
be interpreted as the terms describing the phenomenon 
of multiple reflections between the first and second bar- 
rier. 

The meaning of the picture just outlined is then, quite 
trivial. Ri and T 2 corresponding to the second barrier 
are simply given by the one-barrier coefficients Rob, Tob 
times a phase factor which takes into account the fact 
that this barrier starts at x = L, while the reference 
point in our discussion is at x = 0. Instead, R\ and T\ 
corresponding to the first barrier are given by the sum of 
two terms: the first one is the no multiple reflection co- 
efficients while the second one describes the phenomenon 
of multiple reflections. However, it is remarkable that no 
multi-pie reflection coefficients R\ and T° in Eqs.(29) do 
not coincide with the one-barrier coefficients Rob and 
Tqb- This is an obvious consequence of the fact that 
the scattering probability from the first barrier, neglect- 
ing multiple reflections, cannot be equal to the unity and 
the extra terms Rq and Tq in Eqs. (29) must be present 
in order to achieve the probability constraint in Eq.(3I). 
On the other hand, the scattering probability, including 
multiple reflections, must be equal to 1 (according to Eq. 
(30)), so that we can deduce that the quantities Rq and 
Tq are related to the multiple reflection coefficients i?R 
and Tr. It is very easy to obtain from Eqs. (28) that 
Rq and Tq differ from i?R and Tr just for a phase factor 
(depending on L and a): 



R q 
R^ 



Tr 



_ e ik(L-2a) 



(33) 



Then, multibarricr tunneling is a highly non-local phe- 
nomenon driven by multiple reflections, whose influence 
on the determination of the reflection and transmission 
coefficients is (indirectly) present even in the case in 
which they are neglected. 



ered and introduce the quantities: 

0! ^argjT! e lka } 
02 = arg{T 2 e tka } 
( />s = arg{<? e ^ L - a >} 



(34) 



whose derivatives with respect to frequency give the 
phase times for the first barrier crossing, the second bar- 
rier crossing and the time associated to multiple reflec- 
tions, respectively. Since T = TiT 2 S from Eqs. (23), the 
total tunneling phase is given by: 



= <f>i + {<t>2 - kL) + 4> s 



(35) 



This relation leads to the obvious conclusion that the tun- 
neling time t is the sum of the partial times T\ and r 2 ** 
spent to travel across the first and second barrier, respec- 
tively, plus the time ts required by multiple reflections 
in the inter-barrier region of length L — a. However, it is 
interesting to evaluate the explicit expressions for these 
times and, from Eqs. (34) we get: 



t>o ~ — I- ka 



kL 



(36) 
(37) 



where 4>o = arg{ToBe ifca } is the one-barrier tunneling 
phase time and ^ 



(38) 



While the time to cross the second barrier equals exactly 
the one-barrier tunneling phase time (see Eq.(37)), from 
Eqs. (36) and (38) we see that: 



ii + 0s = , 



(39) 



that is the time spent to travel from the starting edge of 
the first barrier to the starting edge of the second one is 
zero. Something similar to this statement has already 
been suggested in literature (see, for istance, Rcf. [If]), 
but now we have a quantitative proof for that. Moreover, 
we can also deduce that, due to multiple reflections, the 



C. Tunneling phase time 

Let us now consider the tunneling phase time r in Eq.(10) 
corresponding to the double barrier crossing here consid- 



|£0|2 + |T 0|2 = l+T 2 e -2 X a ^ 

which makes impossible to give a direct physical meaning to 
R°,T°. 



"Note that the time T2 corresponds to the phase 02 — kl, 
since the travelling along the distance L is already taken into 
account in <f>i + 0s or, in other words, in the expression for 
the coefficient T2 in (25) we have already considered the shift 
from x — to x — L. 

^For opaque barriers, the leading term in S is, from 
Eqs. (22), (25), (26): 

= (fc-Hx) 2 c -ikL/2 2xfc 

4i\k 2xkcoskL/2 - (k 2 — \ 2 ) sinfcL/2 
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time to cross the first barrier is usually different from the 
one-barrier tunneling phase time since: 



(40) 



where 



b Q = avg{T Q } = 20o 
6 R = arg{T fl } = 20 o 



kL — 2ka 



are the phase times corresponding to the terms Tq and 
Tr, the equality holding true only in the case in which 
the inter-barrier distance coincide with the barrier width, 
i.e. L = 2a. 



IV. CONCLUSIONS 

In this paper we have studied the tunneling of a particle 
or a photonic wavepacket through an arbitrary number 
TV of finite rectangular opaque barriers and obtained an 
analytic expression for the total transmission coefficient 
Eq.(9). From this we have confirmed and generalized to 
the present case what was found earlier for a system of 
one [8] or two [11] barriers: the (total) tunneling phase 
time is independent both of the barrier width and of 
inter-barrier distance. These features have been observed 
experimentally for single [2]- [6] and double barrier [10] 
tunneling using photonic setups. 

Amazingly enough, we have further found that, although 
the tunneling probability decreases exponentially with 
the barrier thickness and with the number of barriers (in 
the opaque barrier limit), the tunneling time does not 
depend even on the number of barriers crossed, i.e. it is 
the same for one, two or more barriers. Moreover, when 
considering resonant tunneling, we have also shown that 
the position in frequency (or energy) of the resonance of 
the structure is independent of the number of barriers as 
well. These two novel predictions can be experimentally 
tested using, again, photonic devices. 
In order to obtain a physical picture of what happens 
in the system considered and, especially, of the peculiar 
features outlined above, we have studied the role of multi- 
ple reflections between the barriers on the tunneling and 
found this to be a highly non-local phenomenon. In fact, 
as shown in Sect. Ill, even in the case of increasingly 
large separation between the barriers, the effect of mul- 
tiple reflections cannot be avoided at all. In particular 
multiple reflections play a crucial role in the understand- 
ing of the intriguing results on the (total) tunneling time 
quoted above. Though in Sect. Ill we have dealt with a 
two-barrier system, the main results achieved can be eas- 
ily generalized to multibarrier tunneling as follows. For 
N barriers the partial reflection and transmission coeffi- 
cients corresponding to the first N — l barriers are clearly 
influenced by multiple reflections occurring in the inter- 
barrier regions, while those associated to the last barrier 



are not and coincide with one-barrier coefficients up to a 
phase factor. In particular, as shown in Sect. IIIC, the 
tunneling phase time for crossing only the last barrier 
equals that for a single barrier structure. Since the total 
tunneling time for crossing all the barriers coincides as 
well with the one-barrier time (see Sect. II), we immedi- 
ately deduce that the time for travelling from the starting 
edge of the first barrier to the starting edge of the last 
one is zero. Note that such a result can be achieved only 
if we taken into account multiple reflections and, in any 
case, the partial times for crossing single barriers are usu- 
ally different from the one-barrier tunneling time. 
Finally we point out that our findings also agree with 
the recent results reported in [12], according to which 
a wavepacket travels in zero time a region with TV 
5— function barriers. In fact, as said above, the total 
tunneling time coincide with the transit time for the last 
barrier or one-barrier phase time. From [9] (see Eq.(13) 
of that paper) we then see that, for \ — * 00 , this time 
tends to zero, thus recovering the result of Ref. [12]. It 
would then be nice, in the future, to make the connection 
between multiple reflections studied here and the tunnel- 
ing interpretation in terms of superoscillations quoted in 
[12]. 
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APPENDIX A: COEFFICIENTS FOR N = 2 

From Eqs.(7) we obtain the following expressions for the 
coefficients describing the propagation through two suc- 
cessive opaque barriers: 



R ~ i — — [1 + 2i sin k(L - a) T e~ 2xa ] 



B 



k + i\ 

A 2 ~ (fc ~ 2 smk(L - a) T e^ a 

k — i\ 2\k 

k — %x \ 2k [ 1 (k — ix) 2 



2k 



k + ix [k-ix 
■ T e- 2 * a ] } 



1 



2 x k 



sinfc(i — a) 



A 3 ~ e- lkL T e~ xa 



k + ix 
At ~ 

2k 

B 4 ~ 
T ~ 



k l X n ikL j: e - X a 



k + ix 
4ixk 
(k + ix) 



T e~ xa 
- T e~ 2xa 



(Al) 
(A2) 

)• 

(A3) 
(A4) 

(A5) 

(A6) 

(A7) 

(A8) 
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(in all these expressions we have neglected terms of third 
order in e xa ). 
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